Temperature Dependence of Even-Odd Effects in Small Superconducting Systems 
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Generalized BCS equations which consistently include the projection on the particle-number parity 
are derived from a systematic variational method. Numerical solutions are given that are illustrative 
of ultra-small metallic grains. Compared to the BCS approximation, the pairing is slightly enhanced 
in even and noticeably inhibited in odd systems. In the latter case, the gap increases with T near 
T = 0; for a suitably small pairing this can result in a reentrance phenomenon. Paramagnetic 
effects are evaluated as functions of the temperature and the magnetic field. 
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The qualitative differences between the properties of 
even and odd nuclei have long been known, as well as 
their explanation by pairing correlations Even-odd 
effects have also been displayed in recent years in super- 
conducting metallic islands |Q and, during the last two 
years, in ultra-small Aluminium grains In the nu- 
clear case, the ground and lowest excited states are suc- 
cessfully described by the BCS theory and its Hartree- 
Fock-Bogoliubov (HFB) generalization Indeed, at 
zero temperature, the BCS state is a superposition of 
configurations with even particle number only, which al- 
lows the currently used models to discriminate between 
odd and even systems. However, in the case of hot nuclei 
(which will be investigated intensively by means of the 
next generation of 7-detector arrays) or of small metal- 
lic superconducting systems, the non-zero temperature 
BCS Ansatz mixes configurations with even and odd 
particle numbers that have dissimilar physical proper- 
ties. This additional spurious mixing requires projecting 
the thermal BCS state onto subspaces having, at least, a 
well-defined number-parity (if not a well-defined particle 
number). Several theoretical studies have already been 
devoted to this question ||]. 

We rely here on a general variational scheme to 
derive an extension of the BCS theory which consis- 
tently incorporates the number-parity (A^-parity) projec- 
tion. The method is devised to optimize the (unnormal- 
ized) characteristic function i^(^) = lnTrP^e~'^^A(^), 
with K — H — jiN and A(^) = exp(— ^^Q-y) where 

H, N, and ^-y are the Hamiltonian, the particle- 
number operator, the observables of interest and the as- 
sociated sources. The trace (Tr) is taken over the full 
Fock space. In the unnormalized equilibrium density op- 
erator P^e~^^ , the factor is the projection 



Pr. 



on the subspace with even {rj = +1) or odd (77 = — 1) 
particle number. To determine variationally (^(0 we in- 
troduce H the action-like functional 



Z{p(r),i(r)} =TrP(/3)exp(-^e7Qi 



(1) 



where the trial quantities are the operators A{t) and 
27(t). The density-like operator P(r) is meant to include 
the projection factor P,, and to satisfy the initial bound- 
ary condition P(0) = P,,. For unrestricted trial spaces, 
the stationarity conditions of the functional (P yield the 
(symmetrized) Bloch equation for P(r) and, redundantly, 
the Bloch analogue of the backward Heisenberg equation 
for A(t) with the final boundary condition A{P) = A{^). 
When the trial spaces for T>{t) and A{t) are restricted, 
the two equations in general become coupled and the sta- 
tionary value of the functional (^ provides an optimum 
/orexpv5(^). 

We are only concerned in this letter by the evalua- 
tion of (projected) thermodynamic quantities. These are 
given by the limit = (i.e., A{0) — 1); indeed (f{0) 
is the thermodynamic potential associated with the P^- 
grand canonical equilibrium in which the particle number 
can fluctuate but not the A'^-parity. In this limit and for 
the trial spaces (||) that we shall use, one can check that 
the optimization of is equivalent to the usual mini- 
mization of the thermodynamic potential . 

For the Hamiltonian H we take, as in ||^, the schematic 
BCS form 

= ^ ep{alap + a^ap) - P M - ^ Gpga^alaqag , 



in which aj, and Op are the creation operators associ- 
ated with the paired single-particle states p and p. The 
magnetic-moment operator M = J^p ~ 0'j,(^p)i cou- 

pled with the magnetic field B (in reduced units), gives 
rise to paramagnetic effects; for rotating deformed nuclei, 
the term BM would be replaced by loJz where lo is the 
angular velocity and the z-component of the angular 
momentum. The pairing coefficients Gpq are real, sym- 
metric and, for simplicity, the diagonal elements Gpp are 
assumed to be zero. 
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As variational Ansatz we take 



P(t) = P^e- 



-tW 



Mr) 



-(l3-r)W 



(2) 



with, for the trial operator W , the standard form 



-r ep6]j6p) , (3) 



where 6j and bp are the usual quasi-particle operators 
(hp = UpQp + VpOp, bp = Upttp — Vpaj^). In the static 
problem we are considering here, the quantities Up and 
Vp can be taken real and non- negative, with Up + Vp = 1. 
Thus our independent variational parameters are the real 
numbers Cp, ep, Vp and ho- We get an explicit expression 
for the functional (|l|) by inserting Eqs. (Qj^) and using 
the generalized Wick theorem for the two terms gen- 
erated by the projection i^,,. New projected BCS equa- 
tions are obtained by writing the stationarity conditions 
for this expression. The familiar BCS equation for the 
gap is replaced by 



Ap - - ^ Gpq^{tq + tq)- 
9 ^ 



1 + r]ro{tptp)-^ 



(4) 



where the quantities tp, tp, tq, Cg, are defined by 
tp = tanh(i/3ep) , tp = tanh(i/3ep) , '^o = J]^ tptp 



Up"^ + Vp^ = 1 and by the relation 2upVp / (up"^ — Vp^) = 
lS.pl [ep — /i). All the quantities Ap, Cp, Cp, Cp, rg = 
Tr e^'^^e-'^^/Tr e"''^, Xp, fci, Up, Vp depend on the A^- 
parity index rj = ±1. The standard BCS equations are 
recovered for r/ — 0. In Eq. (||) the projection effects 
occur indirectly via Cp and ep, and directly through the 
last fraction. The form of this fraction entails : (i) a pre- 
dominant influence of the states near the Fermi level, (ii) 
odd (even) projected gaps Ap smaller (larger) than their 
BCS counterparts, and (iii) greater deviations from BCS 
for the odd than for the even systems. 

We further assume, as is often done, that the twofold 
degenerate single-particle energies ep are equally spaced 
with a level density wy- We use a cut-off energy A/2, 
above and below the Fermi level, such that wfA = 100. 
The pairing coefRcients Gpq are chosen of the form 



G(l 



The projected and the BCS gaps depend 



on p ; however, the variations of Ap with p never exceed 
a few percent in our calculations. A convenient energy 
scale is provided by the BCS gap A associated, for T — Q 
and 5 = 0, with the quasi-particle of lowest energy. This 
model, depending on the values of wpA, or alternatively 
of wfG, schematically describes superconducting islands, 
heavy nuclei or ultra-small metallic grains. The corre- 
sponding variational equations above have been solved 
in Ref. |^ for B — 0. The results display the growing 
magnitude of the even-odd effects as wpA decreases. In 
view of recent experimental advances [^j , we give prefer- 
ence here to parameters typical of ultra-small grains. 



(ep - m)2 + Ap 



(5) 



The quasi-particle energies Cp and ep appearing in the 
thermal factors tp and tp differ from the BCS expression 



indeed, ep satisfies 



B 



■qr^t 



tp / Xp{tp -l- ip ) 



■(i 



tp'' - Tyro V 1 + Waihh) ^ 
where Xp and ki are given by 



^ + 2fci) , (6) 



A„A 



pq , ^ 

e«e 



Wo + Wo 



pcq 



(tq+tq-tq-tq) 



l-tp^ 
2 I 



"4 $Z ^pff^^P 



+ P < 



P) 



Wo 



p ^ p 



<1 



<1 




The quasi-particle energy ep is obtained from Eq.(g) by 
the exchanges tp tp and B <-» —B. Altogether, one 
has to solve numerically the self-consistent Eqs. (^) 
and (|^) for Ap, Cp and ep. The coefficients Up and Vp 
are then given explicitly by the normalization condition 



FIG. 1. Even/odd projected gaps. The gap Ap is plotted 
(for B = 0) as function of T for various values of topA, or 
equivalently of wfG- The upper part (a) corresponds to the 
even system (N) = 100 and the lower part (b) to the odd 
system (N) = 101. 
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For B = and for various values of wpA, Fig.l shows 
the variation with T/Tc (Tc denotes the BCS critical tem- 
perature for i? = 0) of the ratio Ap/A, where Ap is the 
temperature-dependent (odd or even) gap Ap associated 
with the quasi-particle of lowest energy Cp. According to 
Eq. (Q), both the gap values and the critical temperatures 
show that pairing is enhanced (inhibited) by projection 
on even (odd) particle number. In the odd case, near 
T = 0, a new phenomenon appears: the growth of the 
pairing correlations with temperature. This effect can be 
understood as follows. At T = the unpaired particle 
occupies one of the two states eo — — /i, precluding 
them to participate to the pair components. For (small) 
increasing values of T, thermal excitation induces the 
migration of this unpaired particle to higher states, thus 
making the strategic states at the Fermi level available 
to the formation of pairs. 
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FIG. 2. Dependence upon T of the gap Ap for wpA = 1.19 
{wfG = 0.233) and for magnetic-field values such that 
wfB=0- 0.2; 0.4; 0.6; 0.8. The three plots correspond to 
(a) the even-iV projection {{N) = 100), (b) the BCS theory 
{(N) = 100), and (c) the odd- projection {{N) = 101). Note 
the scale difference between (c) and (a). 

For suitable values of wpA (in our example 1. < 
wpA < 1.14), this increase of pairing correlations ver- 
sus T (near T — 0) can even generate a reentrance phe- 
nomenon with two transition temperatures (Fig. 1-b). 
In such a case pairing is absent at T = ; as T in- 
creases, it switches on at some low critical temperature, 
then switches off at the higher familiar critical temper- 
ature The two critical temperatures merge for 
wpA « 1, which corresponds to wpG « 0.225; for weaker 
values there exists no superconductivity in odd systems, 
whatever the temperature. In contrast, in even systems. 



pairing does not disappear until wpG < 0.182. 

The comparison of (a) and (b) in Fig. 2 shows that, at 
all temperatures, the even-A^ gap resists magnetic fields 
better than the BCS gap. It is also seen on Fig. 2 that, at 
T — Q, the gaps Ap do not depend on the magnetic fields; 
this is because the field intensities are such that B < A. 
When B > A, for the value wpA « 1.19 of Fig. 2, pair 
correlations disappear completely. By comparing (b) and 
(c) one sees that the odd-A^ gap is significantly smaller 
than its BCS counterpart and more easily distroyed by 
the paramagnetic effects. Nevertheless, the growth with 
T of the odd-A^ gap at small T persists in the presence 
of a sufficiently small magnetic field. 

Once the self-consistent equations (^ and (^ are 
solved, the projected average particle number, magneti- 
zation, energy, entropy (and hence the other thermody- 
namic quantities) can be obtained from the formulae 



{N) 

m 



1 



2(1 + Tyro) 



s 



BCS 



(^tp -f- tp) 



2(l + ^ro) 

Ap" 
X 

1 -t- riroitptp''^^ 



- B (M) 



2(l + ryro) 



2(1 + Tyro) 

E[(^p^ -ip)ep+P^p] 



+ln[i(l + ,7ro)] . 

These quantities satisfy the fundamental thermodynamic 
relations, as warranted on general grounds by the varia- 
tional nature of our approximation. Various limits (low 
and high temperatures, continuous level density, large 
volume) and numerical illustrations are worked out in Ig] 
for B = 0, showing the qualitative differences between 
even and odd systems in physical quantities such as the 
specific heat. 

Here we show on Fig. 3, for the same values of w-pA 
and of the magnetic fields as in Fig. 2, the magnetization 
(M) versus the ratio T/T^ for the BCS theory, for the 
even (A^ = 100) and for the odd (A^ = 101) projected 
cases. All the curves converge to the Pauli paramagnetic 
values 2wfB at high enough temperatures. Until the in- 
dependent transition temperature is reached, the BCS 
magnetization is always larger than the even-A^ one, but 
both behave similarly. The behavior at small T of the 
odd-A^ curves is markedly different. While for B — 
the magnetization is zero as expected, any value of B, 
however small, lifts the single-particle degeneracy of the 
spectrum and causes the zero-temperature value of (M) 
to take integer values. For the small magnetic fields con- 
sidered here {wpB < 1), one has (Af) = 1 for T = 0. 
As T grows, (M) increases or decreases in odd systems, 
depending on the value of B, so as to reach eventually 
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the Pauli paramagnetic value. 



A 




FIG. 3. The magnetization (M) versus T /Tc for dif- 
ferent values of the magnetic field B and for uipA = 1.19 
(tiJpG = 0.233). The dashed lines correspond to the even-A'' 
projection ({TV') = 100), the solid lines to BCS {{N) = 100) 
and the dotted lines to the odd-A^ projection ((iV) = 101). 

Our curves exhibit a spurious non-analytic behavior at 
the critical points. This is a consequence of the approx- 
imation (|) which, hke the BCS one, gives rise to either 
a broken or an unbroken invariance with respect to iV, 
depending on the temperature. In the exact treatment 
of a finite system without sources, the order parameter 
should always be zero and there can be no discontinu- 
ity at any temperature. However, allowing an order pa- 
rameter while restoring the broken symmetry may re- 
produce the rounding-off of the transition For in- 
stance, in the present problem the difficulty above may 
be removed by replacing by the projection Patq = 

(27r)~^ d^ exp(i6'[iV — A'o]) on the particle number 
iVo, thus restoring the invariance broken by W . However, 
in this case the trial Ansatz V{t) — Pis!qT{t)Pis!q, where 
T(t) is the exponential of the most general quadratic 
form in a J and aj, must include projections on both sides 
1^. (In contrast, the Ansatz (^ contains a single projec- 
tion since Pr^ commutes with W .) After the substitution 
of Prj by Patq, which amounts to replacing a two-term 
sum by a double integral over 0, the functional (|^) re- 
mains an operational tool since one can still use Wick's 
theorem. The optimization of (|l|) then provides a gen- 
eralization of the HFB theory which takes properly into 
account the particle number. In this scheme, where vari- 
ation is performed after projection, the quantities Up Vp 
are expected to be non-zero at any temperature, and to 
yield smooth transitions. Likewise, in (cold) deformed 
nuclei the evolution of the moment of inertia with re- 
spect to the angular momentum is smoothened by the 
P/Vo-projection, in agreement with the data |l^ ]. 

Projected fluctuations and correlations could be eval- 
uated, along the lines of Sect. 4 of through the ex- 
pansion of ip{£,), as given by (|^), up to second order in 



the sources ^-y . Finally, as proposed in |14| , transport 
properties could be worked out by combining the static 
variational principle (|l|) with a time-dependent one. 
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